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ABSTRACT 

We present a new study of Ha/FUV flux ratios of star forming regions within a 
sample of nearby spiral galaxies. We search for evidence of the existence of a cluster 
mass dependent truncation in the underlying stellar initial mass function (IMF). 
We use an automated approach to identification of extended objects based on the 
SExtractor algorithm to catalogue resolved H ii regions within a set of nearby spiral 
galaxies. Corrections due to dust attenuation effects are applied to avoid artificially 
boosted Hg/FUV values. We use the BPASS stellar population synthesis code of 
lEldridge fc Stanwavl (|2009D to create a benchmark population of star forming regions 
to act as a reference for our observations. Based on those models, we identify a zone 
of parameter space populated by regions that cannot be obtained with a cluster 
mass dependent truncation in the stellar IMF imposed. We find that the investigated 
galaxies display small subpopulations of star forming regions falling within our zone 
of interest, which appears to be inconsistent with imposing an IMF truncation at 
maximum stellar mass dependent on the total cluster mass. The presence of those 
regions is expected in models using both stochastic and sorted sampling of the full 
extent of the stellar IMF. This result persists after taking dust attenuation effects 
into account. We highlight the significance of stochastic effects in environments with 
low star formation activity and in studies describing systems associated with small 
physical scales. 
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1 INTRODUCTION 

The stellar Initial Mass Function (IMF) is the key prop- 
erty characterising a newly formed stellar population. It de- 
scribes the underlying probability distribution of the forma- 
tion of a star within a particular mass range. Currently, the 
most commonly used parametri sation is a two -part power 
law £,{m) oc m~°'' presented by iKroupal (|200ll ). with ai = 
1.3 for 0.1 ^ m/M(:y ^ 0.5 an d a2 = 2.35, identical with 
the slope from ISalpeteiJ (|l955l ). for 0.5 ^ m/MQ ^ n^maxi 
where m^^x is the fundamental upper stellar mass limit. 
At Solar metallicity m^ ax = 150Mq is widely adopted 
(jZinnecker fc Yorkd 120071 , and the references within) , al- 
though 120Mq is also often used. 

Over the recent years a considerable amount of dis- 
cussion was devoted to the question of whether the stel- 
lar I MF is indeed u niversal as postulated by many w orks 
(e.g. iKroupalliooH) . iBastian. Covey, fc Meveil (|2010l ) and 
the references therein provide an extensive review of the 
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observational evidence on IMF variation. The Integrated 
Galactic IMF (IGIMF) t heory l|Kroupa fc WeidneiJlJOOsI : 
IWeidner fc Kroupal l2005l. 20061) has been de v eloped to ex- 
plain the claims (e.g IPrisinzano et al.l l2003l : iKalirai et al.l 
i2003l : iMoraux et al.ll2004 ) for a presence of a top-Hght IMF 
(i.e. declining more steeply within the high mass regime). 
The core of the IGIMF paradigm is that the effective galac- 
tic IMF is composed of the contributions from individual 
star forming events. Each of those events corresponds to 
formation of a cluster and a separate sampling of the ade- 
quate underlying stellar IMF. Note that this view encom- 
passes both clusters that remain bound and that disperse 
to form stellar associations. 

The semi-analytic implementation of the IGIMF the- 
ory is founded upon two underlying assumptions. First, the 
maximum mass of a cluster formed in a galaxy (Mgci.max) 
is taken to be a functio n of the star fo rmation rate (SFR) 
of the galaxy ( Wcidnor . Kroupa. fc Lar sen 2004). Second, 
the maximum mass of the star rrimax that can be formed 
within a newly formed cluster of mass Med is dependent 
on Mcci itself. This semi-analytical relation was presented 
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in IWeidner fc Kroupal ||2006|) and was described i n furt her 
detail bv lPflamm-Altenburg. Weidner. fc Kroupal (|2007l , af- 
terwards PA07). We will be referring to it as the PA07 re- 
lation. 

Two compilations of published Milky Way clus- 
ters datasets attempte d to d emonstrate this relationship. 
iMa schbcrgcr & Clarkg ()2008l ) analysed a sample of Milky 
Way clusters and found no significant departures from the 
stochastic sampling of the full extent of the IMF in low 
mass clusters, as described earlier by lElmeereenl lJ2006f ). 
No alternative sampling methods have been investigated, 
though. Weidner, Kroupa, & BonncU (2010) presented a 
similar compilation of a larger cluster sample extending to 
higher masses, containing observational evidence support- 
ing the PA07 relation. Those authors concluded that a re- 
lationship between Mod and mmax is a physical constraint 
rather than a statistical effect. They also point out that in- 
vestigating a wider cluster mass range is necess ary for dis- 
tingui shing between different sampling methods. lElmeereenl 
l|2006r ) interprets the relationship between Mgci and mmax 
as a statistical connection originating from the size of sam- 
ple effects. 

Rather than attempting to verify the PA07 re- 
lation directly, one can search for its manifestations 
on galactic and sub-galactic scales. The semi-analytic 
IGIMF paradigm predicts a divergence of Ha and far 
ultraviolet (FUV) fluxes of galaxies with low SFRs 
(jPflamm-Altenburg. Weidner. fc Kroupal I2OO9I '). This is be- 
cause below a certain SFR threshold formation of stars mas- 
sive eno ugh to have significant ionizing flux becomes im- 
possible. IPflamm-Altenburg et al.l (|2009f ) predict that this 
should translate into a discrepancy between the star for- 
mation rate s recovered from Ha and FUV calibrations. An 
attempt bv lLee et al.l ()200g l to observe this trend for inte- 
grated galactic fluxes has confirmed a systematic turnover 
in Ha/FUV ratio at low SFRs, but also showed that the 
origin of the trend was inconclusive. 

The Afcci - mmax relationship presented by PA07 has 
been widely overinterpreted as an imposition of a strict 
cluster mass dependent truncation in the stellar IMF. 
IPflamm-Altenburg et al.l (2009) predictions have sparked 
further investigations into the potential presence of a 
sharp IMF cutoff ?nmax < nt^a x at lo w clu ster masses. 
iFumagafli. da Silva. fc KrumholJ l|201lh and IWeisz et all 
l|2012h demonstrated discrepancies between the observa- 
tions and predictions from a variable IMF t runcation com- 
bined with the Meci.m ax - SFR relationship (jWeidner et al.l 
I2OOJ) . lEldridgd (|2012l ) has found similar discrepancies when 
using a truncated IMF in absence of a Mcci,max-SFR depen- 
dence. 

Previous studies of integrated galactic properties ques- 
tion the validity of using a truncated IMF models inspired 
by the semi-analytic IGIMF. However, they cannot deter- 
mine whether the inconsistencies arise due to the implemen- 
tation of the IMF truncation, the Mod, max - SFR relation- 
ship or both. The aim of this study is to verify the first 
of those options by investigating star formation on sub- 
galactic scales. To accomplish this, we perform a search 
within a sample of nearby spiral galaxies for a popula- 
tion of star forming regions that are inconsistent with the 
presence of a cluster-mass-dependent upper truncation to 
the IMF. We use the Binary Population And Spec- 



TRAL Synthesis (BPASS) stellar population synthesis code 
JEldridge fc Stanwavll2009l) to extend the Hq/FUV flux di- 
vergence predictions of IPflamm-Altenburg et al.l ([20091) to 
the regime of individual clusters with various approaches 
to IMF sampl i ng. T his study is complemen t ary t o work of 
ICorbelh et ^(l2009h o i iM33 , ICalzetti et al.l(|201ol l on NGC 
5194 and iKoda et al.l (|2012l ) on the extended UV (XUV) 
disc of NGC 5236. 



2 DATA SAMPLE 

For the purposes of this experiment we have selected a sam- 
ple of seven nearby nearly face-on spiral galaxies for which 
ground-based Ha, FUV, NUV and Spitzer Space Telescope 
MIPS 24/im data are available. The galaxies were selected to 
span a narrow distance range. This is because more distant 
galaxies are not resolved enough to infer the properties of 
small enough star forming regions. On the other hand, the 
assumption of spatial coincidence of the sources of ionizing 
radiation and Ha emission breaks down in better resolved 
nearby galaxies 

Table [l] provides an overview of the properties of 
these galaxies. The galaxies' positions, morphological types 
and E(B — V) values were obtained directly from the 
NASA/IPAC Extragalactic Datab ase (NED) . In mo st cases 
we adopt the same distances as iLee et al.l (|201ll ) did in 
their study. The one exception is NGC 628, where we use 
more recent measurements based on the observed luminos- 
ity functions of planetary nebulae. Throughout the analysis 
we adopt the distance uncertainties as quoted in the rele- 
vant papers. The references are listed in the caption of Table 

m 

The MIPS 24^m data are taken from two multi- 
wavelength data sets: Local Volume Legacy survey (LVL, 
iDale et al.ll20p9|) and Spitzer Infrared Nearby Galaxies Sur- 
vey (SINGS. iKennicutt et"a l. 2003). Most of the narrow- 
band Ha images are obtained from sets of ancillary ob- 
servations to the two abovementioned surveys. In the case 
of NGC 300 we used data from t he Survey for Ioniz ation 
in Neutral Gas Galaxies (SINGG. 'Meure r et al.|[200^ ). For 
each of the galaxies we selected the highest quality Ha data 
available. The references for Ha images are quoted in the 
description of Table [l] Note that due to the large mismatch 
between the point spread functions (PSFs) of FUV and in- 
frared (IR) data, high angular resolution of the Ha data is 
not necessarily a significant a dvantage. 

We used GAIA software IjJenness et al.ll2009r ) to manu- 
ally remove star subtraction artefacts from the narrow-band 
Ha images using the inbuilt patching functions. We typi- 
cally used circular apertures with a third order fit (noise es- 
timate scaling factor of 1.0) to a background annulus in the 
patching procedure. The annuli's radii were selected in or- 
der to avoid neighbouring objects. GAIA fits a single back- 
ground surface across all the regions selected. Hence, in or- 
der to optimize the quality of the process, patching was per- 
formed for each region individually. In more complex cases, 
such as crowded regions or presence of artefacts associated 
with prominent diffraction spikes in the continuum image, 
appropriate patching region shapes were selected to com- 
pensate. In such situations, we were often forced to abandon 
the use of background annuli in favour of individually as- 
signed sets of background regions. We find that as a result 
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Name RA Dec distance typo E{B — V) log(LHQ) log(ipuv) telescope 

(J2000) (J2000) [Mpc] [mag] [erg s^i] [erg s^i] 

NGC 300 00:54:53 -37:41:04 2.0(i' SA(s)d 0.0130 40.18 42.46 Danish 1.54 m(5) 

NGC 628 01:36:42 +15:47:01 8.6(2) SA(s)c 0.07 40.87 43.02 CTIO l.SmW 

NGC 2403 07:36:51 +65:36:09 3.22(1' SAB(s)cd 0.04 40.78 42.83 KPNO 2.1m('5) 

NGC 3031 09:55:33 +69:03:55 3.63(1) SA(s)ab 0.08 40.77 42.76 KPNO SchmidtC^) 

NGC 5236 13:37:01 -29:51:56 4.6l(^) SAB(s)c 0.066 41.25 43.21 CTIO 1.5m(*) 

NGC 5457 14:03:13 +54:20:56 6.7(1' SAB(rs)cd 0.01 41.33 43.61 KPNO SchmidtC^' 

NGC 7793 23:57:50 -32:35:27 3.9l(*) SA(s)d 0.02 40.58 42.67 CTIO 1.5m('') 



Tabl e 1. Overview o f the galaxies sample. Ho luminosities taken from iKennicutt et al.l 1 120081) . I/puv = -^FUV^FUV values taken 
fr omlLee et al.l 11201 ll). R efer ences indicated in the " tele scope" column are the sources fo r Ha data us ed i n the analysis. Ref e rence s: 



(1 Freedman et alj 1 



■ Kennicutt et al.l 1 



2001) 



20031) 



(^ JHerrmann et alj 1I200I): (3^Sa.ha et al.l 1I2OO6I'): ("^IP ictrzvhskr et alj 1 I2OIOI) : (5) |Larsen fc Richtlej lll999l 'l: 
C^ teoopes. Walterbos. fc BothunI 1I2OOII ): ( ^Meurer et aU (.2006,) : 
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Figure 1. NGC 2403: continuum subtracted Ha emission prior 
to patching out of artefacts. 



of patching, the quality of object identification is signifi- 
cantly improved. We exclude the regions near the patched 
out objects from the output of our search algorithm. 

Where available, we retained the published continuum 
subtractions. This was not the case with NGC 2403 ancil- 
lary data from SINGS, where astrometry of the published 
R-band continuum image was incorrect. We used GAIA 
software to rederive astrometry for both imag es by match- 
ing t he stars to the Guide Star Catalog II rLa sker et al.l 
120081 ). We then subtracted the R band continuum image, 
scaled by a uniform factor ka, off the narrow band Ha im- 
age to correct for the presence of continuum. R-band image 
has been converted to the pixel scale of Ha data and con- 
volved with a 2D Gaussian kernel with FWHM — 3.0 pix (to 
match Ha image PSF) prior to performing the subtraction. 
The value of kn has been derived by performing photometry 
of a sample of 20 manually selected stars on both images 
and adopting the median scaling. However, we note that 
there is a second order effect due to mosaicking that can 
contribute errors of the order of 5 per cent in flux. The fi- 
nal Ha map is displayed in Figure [T] In the case of NGC300 
we us ed the astrometry.net web-based program l|Hogg et al.l 
I2OO8I) to provide an astrom etric calibrations for the images 
of lLarsen fc Richtled l|l999D . 

The FUV images are taken from the LVL survey 
iDale et al.l (|2009l '): [Lee et all l|201ll '). Most of these observa- 



tions were originally presented in the Gala xy Explorer Satel- 
lite (G AL EX) atlas of nearby galaxies bv lCil de Paz et al.l 
112003), which provides the details of the image process- 
ing pipelines used. We cross-checked our flux calib r ations 
against the integrated galaxy values of iLee et al.l ((20111) 
and found excellent agreement. Foreground extinction cor- 
rection was applied adopting the E (B — V) values from 
ISchlegel. Finklpeiner . fc DavisI (119981 ) and the Milky Way 
extinction curve of Pci (JJ392). 



3 METHODOLOGY 

This section describes the methodology adopted in our 
study. First, a summary of our approach to object identifi- 
cation is given. Second, the method of correcting for dust 
attenuation is described. 



3.1 Object identification and photometry 

In or der to identify our regions of interest, we run S Extrac- 
tor (JBertin fc Arnoutslll99(]| ) package using Ha images of 
galaxies as reference to select the regions of recent star for- 
mation. The code first fits a global background profile us- 
ing K.cr-clippin g, similar to that employed in DAOPHOT 
(|Stetsonlll987l ). SExtractor initially identifies objects to 
be groups of pixels with values exceeding a set threshold 
above this background fit. 

Multiple isophotal deblending is applied to the initial 
detections to evaluate significance of the individual peaks. 
The expected contributions from bivariate Gaussian fits 
to deblending objects are then used to allocate contended 
pixels. SExtractor then eliminates spurious detections 
within wings of shallow profile objects by evaluating the 
contribution to the mean surface brightness from the gaus- 
sian extrapolation of neighbouring objects. If the corrected 
surface brightness exceeds the detection threshold, the ob- 
ject is accepted into the catalogue. Typical SExtractor con- 
figuration parameters used in our study are listed in Table 
[2] They were chosen to maximize the detection efficiency of 
the population of faint objects. 

The program is run in single image mode on Ha data 
and in dual image mode on the remaining three bands. Ob- 
jects identified in this way are referred to throughout the 
paper as "Ha selected" . A second set of runs is performed 
using FUV image as the reference, producing what we refer 
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parameter name 



set 1 



set 2 



DETECT.MINAREA 

THRESH.TYPE 

DETECT.THRESH 

ANALYSIS.THRESH 

DEBLEND.NTHRESH 

DEBLEND.MINCONT 

CLEAN 

CLEAN.PARAM 

BACK.SIZE 

BACK.FILTERSIZE 

BACK.TYPE 

BACKPHOTO.TYPE 

BACKPHOTO.THICK 



10 


35 


RELATIVE 


RELATIVE 


5.0 


5.0 


5.0 


5.0 


64 


64 


0.00001 


0.00001 


Y 


Y 


10.0 


10.0 


35 


100 


6 


6 


AUTO 


AUTO 


LOCAL 


LOCAL 


10 


24 



Table 2. Adopted values of SExTRACTOR input parameters. Set 
1 has been used for analysis in native FUV resolution and for 
the analysis in H« resolution in cases of images taken with the 
CTIO 1.5m and KPNO Schmidt telescopes, due to their compar- 
atively poor resolution. Set 2 was adopted for higher resolution 
Ha images. 



to as the "FUV selected" sample. This paper focuses on the 
results from the Ha selected sample. 

Aperture photometry is performed on the list of detec- 
tions by applying the APER procedure of the IDLASTRO 
package. We use the SExtractor defined isophotal aper- 
tures with signal-to-noise threshold of 5.0 and an appropri- 
ate minimum area. We mask out regions that are identified 
outside the galaxies' disks, within their bulges and near the 
artefacts of continuum subtraction that were removed from 
Hq images during the patching stage. A set of circular aper- 
tures is used to allow for spatial offsets between the FUV 
and Hof emission. The aperture selected for any given region 
is the smallest of the set of available apertures that has an 
area larger Wthan the isophotal area fitted to the region. 

We use two sets of ten apertures each. For each gala:xy, 
we choose one of the two aperture sets, the choice being 
dependent on the angular sizes of discernible Hll regions 
due to the host galaxy's distance and telescope resolution. 
Smallest aperture diameter in "set 1" is 4 pixels and the 
sizes increase in 2 pixel steps until the maximum of 24 pix- 
els. "Set 2" apertures range from 10 to 55 pixels with a 
5 pixel increment in diameter. "Set 1" has been used for 
NGC 3031 and NGC 5457, which have low resolution Ha 
data from the KPNO Schmidt telescope. Remaining galax- 
ies had "set 2" applied in the analysis pipeline. 

All the images have been rebinned and realigned to 
the pixel scale of the relevant Ha image using a modified 
frebin IDL routine, with an added fiux conservation op- 
tion. Aperture corrections have been derived by using the 
same photometric approach on rebinned FITS images of the 
PSFs of relevant instruments. The PSF information was ob- 
tained directly from the websites oi GALEX and Spitzer 
missions. ATV package (|Barthll200ll ) was used for quick 
check photometry and result verification, since it allows for 
direct application of the APER routine. 



3.2 Dust attenuation correction 

Dust attenuation produces a strong effect on the locus of 
points on Ha/FUV diagrams. The data points are displaced 



along a vector with a direction fixed by the adopted at- 
tenuation curve and length determined by the magnitude 
of the attenuation. It is universally acclaimed that FUV 
is more significantly attenuated than the Ha nebular line 
(e.g. ICalzetti et al.ll2000n . Therefore, any study which fails 
to correct for this effect will find that their obje cts have ar- 
tificia lly boosted values of the Ha/FUV ratio. iKoda et al.l 
l|2012l ) accounted for this effect using a uniform attenuation 
of Av = 0.1 mag. This approach was possible because the 
XUV disc environment they investigated has low dust con- 
tent. Therefore, spatial variations in attenuation would not 
affect the result of their work significantly. 

In this study, we are dealing with regions within the 
optically bright discs o f spir al galaxies. Adopting the ap- 
proach of iKoda et al.l ()2012l ) would result in shifting the 
locus of all the data points in our figures along the over- 
plotted arrows by a fixed amount. However, a significant 
spatial variation in the attenuati on values within a d isk of 
a spiral galaxy is expected (e.g. iRoussel et al.l 120051 ). The 
distribution of dust te nds to be concentrat ed around the 
sites of star formation (jPrescott et al.l 120071 1. Therefore, a 
uniform attenuation value derived from net galactic fluxes 
will most likely underestimate the real attenuation in our 
objects. 

In view of this issue we turn to estimating individual 
attenuation values for each of the identified regions. We 
adopt a fix ed dust attenu ation curve for all the regions, 
taken from ICalzettil ((20011), yielding the relation Apuv ~ 
1.82 Aho,. 

This scaling relation incorporates both the attenuation 
curve values and the difference between the attenuation of 
stellar light and nebular emission predicted from the typical 
star forming region geometry. The dust corrections are cal- 
culated using the relation between FUV and MI PS 24 nm 
fluxes using the calibration from Section 4.3 of iHao et al.l 
l|2011f ). yielding the relation in Equation U]). Throughout 
the paper we use I/fuv = Afuv^fuv- 



= 2.51ogio(l-|-(3.89. 



-^24/11X1 



)) 



(1) 



We are implicitly assuming that both the nebular emis- 
sion and the UV-bright stellar populations within our aper- 
tures are atte nuated by the same dust distribution. The 
adoption of the lCalzettil (|2001[ ) attenuation relation imposes 
the geometry of clumpy dust shell surrounding a stellar pop- 
ulation mixed with ionized gas. A caveat to keep in mind 
is that the relative geometry of dust, gas and stars in the 
small star-forming regions targeted in our study can assume 
arbitrary configurations leading to a degree of stochastic- 
ity in the relative attenuations of starlight and nebular Ha 
emission. However, Calzetti ot al. ( 2005 ) suggests that this 
model serves as a good approximation for small HII regions. 
The question of selecting appropriate dust attenuation law 
is explored in more depth in Section [5.31 



4 THEORETICAL MODELS 

In order to predict the observational signatures of the IMF 
truncation, one has to generate a set of benchmark model 
populations of clusters. A code capable of populating the 
IMF stochastically and considering discrete stars has to 
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be used in order to provide reliable predictions within the 
regime of low-mass clusters. We note that although we use 
the term "clusters" , we do not model the dynamics of those 
objects. In physical sense,they may either be bound clusters 
or unbound stellar associations. 



4.1 Model description 

In this study we use the Binary Pop ulation And Spec- 
tral Synthesis (BPASS) code by lEldridee fc Stanwavl 
([2009) to provide the benchmark models. BPASS is a stel- 
lar population synthesis code based on a library of 15,000 
stellar evolution models computed using the Cambridge 
STARS code. The mode l s used are descri b ed in detail by 
lEldridge. Izzard. fc Tou^ (|2008l '). lEldridgd (|2012l ) recently 
used the package to investigate the effects of stochasticity 
and interpret the attempts o f observat i onal verification of 
IGIM F presented recently by iLee et all (|2009l ') ; iLamb et al] 
(120101'): ICaizetti et all (|201ol i. We refer the reader to §2 of 
lEldridgd ( 20121 ') for a more detailed description of the meth- 
ods used to create the models, including the implementation 
of binary evolution. The relevant aspects of the model gen- 
eration will be covered briefly here. 

We opted for the use of this package for the follow- 
ing reasons. First, it takes binary evolution effects into ac- 
count while computing the model stellar populations, set- 
ting it apart from the SLUG code used bv iFumagalli et al.l 
(120111'). In light of the recent observational evidence (e.g . 
Pinsonncault & Stanel^ l2006l : iKobulnickv fc FrveJ |20071 : 



Kim inki ct al. 2009) it becomes apparent that binary frac- 
tion among massive young stars is very high. The work of 
ISana et al.l (|2012l ) finds an interacting binary fraction for 
O stars of 69 per cent. Therefore, the effect of interacting 
binaries cannot be neglected if we are after an accurate rep- 
resentation of the properties of young stellar clusters. 

Second, the code produces the spectral energy distribu- 
tions (SEDs) of the populations by integrating over SEDs of 
discrete stars comprising a particular cluster. This approach 
makes the code perfectly suited for an accurate represen- 
tation of the stochastic effects, which become very signifi- 
cant at low cluster masses. BPASS has severa l advantages 
over the models of Bruzual fc CharlotI (12003) which were 
used bv ICalzetti et al.l 1 201(1 ). Already mentioned are the 
stochastic IMF sampling and the inclusion of interacting bi- 
nary stars. One more important difference is the tre a tment 
of Wolf-Rayet (WR) stars. In iBruzual fc CharlotI (|2003l ') 
WR stars are assumed to have atmospheres similar to 
white dwarfs. This incorrectly estimates their contribution 
to the ionizing flux. BPASS uses the detailed grid of cali- 
brated WR a tmospheres, the Potsdam Wolf-R a yet m odels 
presented bv [ Crafener. Koesterke. fc HamannI (|2003) and 



iHamann fc Grafeneij (|2003l . |2004 ). This allows us to have 
greater confidence in the accuracy of estimating the ionizing 
contribution of WR stars in our stellar population models. 
In this analysis we use models at Solar metallicity with 
the metallicity mass fraction of Zq — 0.02. The clusters in- 
cluded in our models span an age range of 1-8 Myr. This 
choice is motivated by the fact that we are expecting the 
objects that would violate the mmax limit to be among the 
youngest clusters within their host galaxies, with ongoing 
star fromation. After consulting published metallicities of 



iMoustakas et al.l (|2010^ we expect Zq to be a good repre- 
sentation for the galaxies in our sample. 

The models cover a mass range from 0.1 to 12OA/0 for 
single stars. This upper cutoff in the model grid is imposed 
by the high mass-loss rates during the main-sequence phase 
at Solar metallicities (|Vink et al.l 1201 ll ). For binary popu- 
lations, the available initial masses are 5 to 120Mq. We 
include every integer mass up to 20Mq , at which point the 
model spacing increases to 5Mq up to 60Mq . Higher mass 
models are computed at 70, 80, 100 and 12OAf0. We note 
that some binary mergers within the models can produce 
stars with masses in excess of 200Mq. Single star models 
have a similar grid of models, with every integer mass up 
to 40Mq. To make an unbiased comparison between the 
single and binary populations we assume the same IMF for 
the primary star masses as the single stars. We then select 
the secondary stars mass from a flat mass ratio distribution 
where, q = ^. Such a mass ratio distribution does not 
affect the overall IMF of the population. The model grid 
in mass ratio have values of g = 0.1, 0.3, 0.5, 0.7 and 0.9. 
In addition we select the initial separation of the binary 
from a flat distribution in log separation. The model grid 
ranges from an initial separation of 10 to lO^'-R© in steps 
of 0.25 dex. Both distributions are si milar to those inferred 
observationally by I Sana et al.l (|2012l ) . Our assumed initial 
distribution leads to approximately two-thirds of our binary 
systems experiencing interactions. For the evolution of the 
secondary in each model we select the most likely result af- 
ter the kick of the first supernova as the evolutionary path 
it will follow. This is estimated using the method of deter- 
mining the binary systems' survival of the first supernova 
event described in lEldridge et al.l (|201ll ). 

Ha fluxes of clusters were derived by calculating the 
number of ionizing photons emitted at wavelengths short- 
wards of 912 A. We assume s tandard Case B recombination 
(|Osterbrock fc Ferlandl 120061 ') so that each lO"'*^ ionizing 
photons give rise to 1 ergs~^ of Ha flux. No flux leakage is 
assumed i.e. all the ionizing photons are absorbed by gas, 
with no corrections for the presence of dust. FUV flux den- 
sities at Apuv were calculated directly from stacking stellar 
spectra in each cluster. 

In order to provide a benchmark for our observations, 
we generated populations of model clusters using three 
methods. The number density contour plots of the result- 
ing models are illustrated in Figure O The left column 
presents models with binary fraction of 100%, the right col- 
umn presents models with no binary systems. In all cases 
we adopt the broken power-law IMF of Kroupa (2001) with 
an absolute maximum initial stellar mass mj^ax = 120Mq. 
The cluster populati ons are generate d in the following man- 
ner, as described in lEldridgd (|2012h . First, a random clus- 



a pc 

iOli 



law cluster IMF w ith a slope of -2 (dc Griis et al.l l2003l : 
iLada fc Ladall2003l ) . The cluster is then randomly filled with 
stars from the Kroupa IMF, truncated at mass-dependent 
JTimax or 771^^^ i whichcver is lower. The stars are added to 
the cluster until the total mass of the stars exceeds the tar- 
get Meci. The last star added is then retained only if the 
mass of the cluster obtained in this way is closer to the value 
of Meci drawn from the cluster IMF. If the total stellar mass 
was closer to target Med without including the final star, it 
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Figure 2. Number density contour plots of BPASS model populations used in this study presented in the Ha/FUV versus Lpuv 
plane. Orange: cluster mass dependent maximum stellar mass (CMDMSM) models corresponding to the truncated IMF; blue: pure 
stochastic sampling (PSS) models; purple: models generated with sorted sampling, preferred by the IGIMF paradigm (see Section |4. II 
for details). Contours are plotted every 0.5 dex in number density. Binning of 0.1 dex adopted in XY. Left-hand side column illustrates 
models with binary evolution effects included. Right-hand side column illustrates models with single stars only. Note that the majority 
of models at low FUV luminosities(< 10^® ergs~^ ) follow the turnoff consistent with the IMF truncation. Points above the dashed line 
are considered to be violating the truncated IMF condition. 



is removed from the cluster. This corresponds to t he co n- 
strained sampling described in lWeidner fc Kroupal (|2006l V 

In the first set of models, we impose an upper mass 
truncation mmax in the cluster-mass dependent form taken 
from PA07. Th is mirrors previous app li cations of th e trun - 
cated IMF by iFumagalli et al.l (|201ll '): iKoda et al.l (|2012l ') 
and others. We refer to this as the cluster mass dependent 
maximum stellar mass (CMDMSM) method. Constrained 
sampling is used. The resulting models are illustrated by 
the top pair of plots in Figure [21 

The second set of models is referred to as pure stochas- 
tic sampling (PSS). We adopt a mass limit TTimax ~ Mod, 



equivalent to allowing star of any mass within our model 
range to form within a cluster, provided it does not exceed 
the mass of the cluster itself. Constrained sampling is used. 
The resulting models are illustrated by the middle pair of 
plots in Figure (2] 

The final set of models is drawn using the same rrimax ~ 
Meci constraint as the PSS population. However, instead of 
constrained samp ling, we apply th e sorted sampling pre- 
ferred by Woidne r fc Kroupal (|2006t ). The difference is that 
upon reaching the target M^ci, the decision between rejec- 
tion or retention is made for the most massive star in the 
cluster rather than the last one to be drawn. The resulting 
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models are illustrated by the bottom pair of plots in Figure 

m 

The synthetic model contours do not cover the param- 
eter space in a continuous manner. This is because we use a 
discreet grid of initial masses and binary parameters. We re- 
frain from interpolating between the grid since non-linearity 
of stellar evolution becomes an even greater issue than nor- 
mally in a binary setup. Instead, for each star drawn from 
the IMF we use the closest mass stellar model. Throughout 
the rest of the paper we use smoothed population contours 
to improve clarity and alleviate the effects of a discrete stel- 
lar model grid. 



4.2 Discriminating between the methods of 
populating the IMF 

The most striking feature in Figure [1 is that CMDMSM 
models display a prominent downturn in the Ha/FUV ratio 
at low Z/Fuv values, mirroring th e galaxy-wide predictions 
of lPflamm-Altenburg et al.l ((20091). As explained before, this 
effect is a result of the lack of massive stars with significant 
output of ionizing flux. 

The same feature is manifested within the two remain- 
ing populations as well. This is because despite the lack of 
an IMF truncation, drawing a massive star from the IMF 
when forming a low-mass cluster is still unlikely. Therefore, 
the majority of model low-mass clusters will have no mas- 
sive stars. Hence, their Ha/FUV ratio will drop in a similar 
manner as in the CMDMSM case. 

Both PSS and sorted sampling realizations display a 
high Ha/FUV ratio branch. Those are the low-mass clus- 
ters (Med < 10^'^Mq) which contain massive (> 20Mq) 
stars with high output of ionizing photons. The locus of 
this feature is similar in both sets of models. This branch 
is populated by a much smaller number of objects than the 
turnover branch. Sorted sampling predicts a lower relative 
population within this area. 

Note that such a bimodality is expected 

(|Cervifio &: Valls-GabaudI 120031 ) in cases where the 
number of stars contributing to an observable is small 
(< 10). The bimodality persists in BPASS models both at 
different metallicities and with binary evolution switched 
off. Removing binaries results in an even sharper distinction 
between the two branches. 

The dashed line in Figure O represents a candidate se- 
lection cut. Any observed regions placing above this line are 
expected to be incompatible with having a truncated stellar 
IMF. Note that there is a handful of CMDMSM models that 
lie above this selection curve as well. These are the models 
harbouring massive stars formed via binary mass transfer. 
They are outnumbered by at least two orders of magnitude 
by the PSS realizations in this regime though. Those objects 
are not present in the single star model realizations. 

In Figure [2] we see that the differences between sin- 
gle star and binary populations are minor. The effect of 
binary evolution is to rejuventate a population, producing 
more massive stars at later times. This can be seen by ex- 
amining the spread of the cluster distribution at the high 
FUV luminosity end. The single-star population display a 
greater spread. This is due to the ageing of cluster models 
causing a drop in Ha ionising flux originating from a dwin- 
dling population of massive stars. The binary systems are 



able to provide the same Ha luminosity for a longer time, 
which is manifested as the a slight boost of the observed 
Ha/FUV flux ratio. The effect of binaries can also be seen 
in the high Ha/FUV ratio arm, where the preferred locus 
for Wolf-Rayet stars is somewhat changed. Furthermore, the 
arm is also more densely populated in the binary evolution 
scenarios. This is especially visible in the CMDMSM case 
highlighted in the previous paragraph. 

To recapitulate, we decide to adopt the high Ha/FUV 
ratio branch of the PSS distribution as the diagnostic signa- 
ture of formation of stars under a stochastic scenario allow- 
ing for sampling of the stellar IMF beyond the cluster mass 
dependent truncation. We note that no such binary test is 
possible between the PSS and sorted sampling models, since 
they both predict the presence of the high Ha/FUV ratio 
branch. 



5 RESULTS 

In this section, we present the resulting descriptions of the 
star forming region populations within the sample of the 
eight investigated galaxies. We compare the resulting pop- 
ulations with the BPASS model predictions. We discuss the 
observed influence of dust attenuations on the behaviour of 
our data. 

Table [5] contains an excrept from the region catalogue 
for NGC 300. The full catalogues are available in a machine 
readable format in the online version of the paper. Regions 
falling outside galaxy disks and near image areas affected 
by artefact patching (Section 13. ip are not listed. 

In all the figures presented in this section of the pa- 
per, the error bars associated with the data points corre- 
spond solely to the propagation of errors associated with 
the photometric measurements and the uncertainty in the 
measurement of each galaxy's distance as quoted by the 
source. Systematics originating from physical misalignment 
of emission peaks in different wavelengths and variation in 
dust attenuation corrections are discussed towards the end 
of this section. 



5.1 Full sample 

Figure [3] presents the behaviour of the Ha selected region 
sample in each of the investigated galaxies. The data points 
have been corrected for the effects on the dust attenuation 
using Equation ([T|. We chose to plot Ha luminosity rather 
than Ha/FUV in order to illustrate more clearly the be- 
haviour of regions with flux limits rather than confident 
flux measurements. 

Note that the diagnostic wing of high Ha/FUV ratio 
models in the PSS population is now represented as the 
broad horizontal band in Ha . The high luminosity trend 
for both models is equivalent to a constant Ha/FUV ra- 
tio illustrated in Figure [1 The turnover in the CMDMSM 
models is less pronounced in this parameter space. Note that 
the attenuation vector is at a shallower inclination than the 
high luminosity trend. 

The green points represent a manually verified subsam- 
ple of the SExTRACTOR detections, referred to later as the 
"clean" sample. It was obtained in the following manner. 
First, all regions were visually checked against the Ha native 
resolution image. We removed from the sample the regions 
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Figure 3. Ha vs FUV luminosity plots of Ha selected region population. Each panel presents the results for one galaxy. Dust attenuation 
corrections have been applied. Black arrows represent a dust attenuation correction vector corresponding to Ajja = Imag. Blue filled 
contours are PSS models, orange contours are CMDMSM models, as in left column of Figure[2] The 50 per cent completeness limits are 
displayed with dashed lines. Light green points represent the manually verified sample of clean regions. Black points are the remaining 
identified objects. 



© 2013 RAS, MNRAS 000, [D-?? 



Ha to FUV ratios in star forming regions of nearby spirals 9 



ID 


RA 


Dec 


dap 


rfap 


Apuv 


log(-^Ha,corr) 


log{LFUV,corr) 


V 

^ corr 


note 




(J2000) 


{J2000) 


[arcsec] 


[pc] . 


[mag] 


[erg s-1] 


[erg s-i] 






2 


00:54:44.00 


-37:35:18.69 


5.89 


57.10 . 


0.41 


37.36 


38.71 


-1.35 


H 


3 


00:54:41.38 


-37:35:39.49 


5.89 


57.10 . 


0.03 


37.02 


38.58 


-1.56 




4 


00:54:50.16 


-37:35:55.14 


5.89 


57.10 . 


0.02 


37.01 


39.05 


-2.04 




5 


00:54:46.50 


-37:36:32.76 


5.89 


57.10 . 


0.21 


37.33 


38.97 


-1.63 


H 


6 


00:54:44.44 


-37:36:40.24 


5.89 


57.10 . 


0.16 


36.99 


38.76 


-1.77 


H 


7 


00:54:54.67 


-37:36:41.15 


5.89 


57.10 . 


0.17 


37.30 


38.75 


-1.46 


H 



Table 3. Catalogue of regions identified in NGC 300. Full catalogues for all the galaxies are available in the online version of the paper. 
Region ID is not continuous since regions rejected based on their location within the image are assigned IDs as well. The aperture 
diameters (ciap) are quoted both in angular a nd physical scales, using the distances listed in Table [T] Apyy is the attenuation in the 
FUV band derived using the lHao et al.l 1 1201 ll ) relationship. The following columns list dust corrected log(LHti) and log(Lpuv) values. 



logd 



-) is the logarithm of the regions' Ha/FUV flux ratio after the dust attenuation. The regions marked with 'H' in 



the final column are in the hand-picked clean sample. Six columns have been omitted in the printed version for compactness. These 
contain Ha, FUV and MIPS 24/im luminosities and the associated photometric errors. 



excessively affected by crowding or deemed to be a part of a 
larger scale structure. The majority of regions eliminated at 
this step were connected with Hll regions displaying shell 
morphology. 

The second step was the removal of regions without a 
clear match in the FUV image or with a match that was sig- 
nificantly spatially offset. This was to eliminate unreliable 
FUV flux estimates, fn cases where several distinctive H ll 
regions were associated with the same large FUV structure, 
they were rejected unless a matching FUV substructure was 
seen. Regions with no associated FUV peak that were lo- 
cated a large distance from any prominent FUV features 
were retained in the sample. 

The population of identified regions follows the theoret- 
ical distribution closely in the high luminosity regime. There 
is a small number of regions present above the main trend. A 
number of those objects have upper limits on FUV luminos- 
ity, which is consistent with their expected low FUV fluxes. 
Regions with Lfuv upper limits at luminosities exceed- 
ing 10"^^ erg s~^ typically have low photon counts and high 
dust attenuation corrections. The clean sample displays an 
even tighter correlation with the predicted trends. The most 
FUV luminous regions have somewhat lower Hq/FUV ra- 
tio values than those predicted by the models. This could 
indicate a presence of older stellar population within the 
aperture. This effect could also be introduced by overcom- 
pensating for dust attenuation effects (see Section 15.31 for 
more detailed discussion). 

We note that adoption of Ha as the reference wave- 
length means we are not sensitive to tracing the turnover in 
the Ha/FUV flux ratio. Indeed, none of the galaxies present 
evidence of this effect within the Ha selected population. As 
expected, the turnover is manifested in the FUV selected 
populations, which are not presented in this paper. 

In the low luminosity end we see a small number of 
regions that begin to deviate from the CMDMSM prediction 
while remaining in agreement with the PSS models. This 
populations is found within NGC 300, NGC 2403, NGC 
3031 and NGC 7793. This is expected, since these galaxies 
are the least distant ones within our sample. Therefore, they 
offer an enhanced sensitivity towards the identiflcation of 
low luminosity objects. 

Figure|4]illustrates an attempt at statistical description 



of this effect. Model and data populations have been divided 
into four luminosity bins and normalized histograms illus- 
trating the distribution of Ha luminosity (Ha luminosity 
functions, effectively) were made. Stacking has been per- 
formed by adding the normalized histograms for each Lpuv 
bin, weighted by the total Ha luminosity of each galaxy as 
listed in Table [T] Each histogram was then renormalised. 
Note that the histograms of data population do not take 
into account any of the regions with flux upper limits. 

Given the behaviour of CMDMSM and PSS model con- 
tours, we do not expect to be able to distinguish between 
the models in the high-mass (i.e. FUV-luminous) regime. In- 
deed, the histograms of model populations for the three bins 
with logjQ(LFuv[erg s~^]) > 38.45 are very similar. The 
models span the same ranges within Ha luminosities and 
the detailed shapes of the normalised distributions are con- 
sistent with each other. The 38.45 < log]^Q(LFuv[erg s~^]) ^ 
39.45 bin offers a minute hint of an offset between the peaks 
of the distributions of PSS and CMDMSM models, but the 
diagnostic significance of this is minuscule. 

This changes as we investigate the lowest FUV lumi- 
nosity bin. The 36.45 < logio(LFUv[erg s~^]) < 38.45 bin 
encompasses the luminosity range spanned by the diagnos- 
tic high Ha/FUV ratio wing. We see a disparity between the 
peaks of the two model distributions. Note also that within 
this bin the PSS models cover an Ha luminosity range ex- 
tending to values higher than those predicted by CMDMSM 
regions by 0.5dex. 

Figure 13] iUustrates that the data population is consis- 
tent with the models within the highest luminosity bins. 
The luminous outliers seen in the highest Lfuv bin can 
be readily explained when looking at pre-binned Figure [3] 
They correspond to the luminous regions within NGC 5457 
that lie on the extension of the constant Ha/FUV ratio 
trend. Our theoretical models simply do not extend to high 
enough masses to account for them. The low luminosity out- 
liers within I/Fuv bins are adressed in a following section. 

A key result is the presence of observed regions 
above the maximum Ha luminosity cutoff predicted by the 
CMDMSM models. The effect is clearly seen in the two 
lower Lfuv bins. Those regions fall within the Ha luminos- 
ity range predicted by PSS models. 

Figure [5] is constructed in a manner analogous to Fig- 



© 2013 RAS, MNRAS 000, [T]-?? 



10 M.T. Hermanowicz, R.C. Kennicutt and J.J. Eldridge 



en 

o 



-5 



z -2 



en 
O 



-4 
35 



ibAb < ' log,o('LFuv) ^ ^8;45' 



38.45 < logio(LFuv) ^ ^9.45 



:Y 



U 




■i9.4'5' <"log,o('LFuv) = '40!45' 




40.45 < log,o(Lpuv) ^ 41.45 




36 37 38 39 

logio(LHa[erg s"^]) 



40 37 38 39 40 

logio(LHa[erg s"']) 



41 



Figure 4. Hq luminosity histograms for objects in a given FUV luminosity bin. Blue dot-dashed lines represent the PSS model 
population. Red dotted lines represent the CMDMSM model population. Black solid lines represent the Ha selected objects stacked 
from the sample of eight galaxies as displayed in Figure |3] All luminosities are expressed in ergs~^ . Regions with flux limits have been 
excluded from the data histogram. Details of stacking arc explained in Section 15.11 



ure|4l by considering only the manually verified clean sub- 
sample. Overall, the clean subsample reduces the number 
of available regions by « 50 per cent. This becomes most 
significant within the lowest Lpuv bin, with the population 
within the cut reduced from 40 to 16 objects. Despite the re- 
duction in size, the clean sample still exceeds the CMDMSM 
predicted Ha luminosity cutoffs within the two lowest Lpuv 
bins. The consistency with both model populations in the 
high I/Fuv regime is also maintained. 



5.2 Effects of region selection 

A major challenge in investigating resolved populations of 
star forming regions is identification of the regions of inter- 
est. The approach described in Section fS.ll has been selected 
to maximize the identification efficiency. However, there are 
a few caveats associated with using automated region selec- 
tion. 

Our analysis is performed on the star forming regions 
within the optically bright discs of the spiral galaxies. In 
this environment, detection completeness becomes a two- 



dimensional function of the position within the image due 
to the effects of crowding. 

We have performed Monte Carlo tests to study the de- 
tection efficiencies within the galaxies. We have inserted 
artificial objects into the image of each galaxy and checked 
their recovery rate by the pipeline. The test objects were 
point sources of various intensities convolved with the 
GALEX FUV PSF. We placed 20 objects at a time for 
300 realisations per flux threshold tested. Ha fluxes inves- 
tigated ranged from 10'^'^° ergs~^ to IQ''^'* ergs"'^ in steps 
of 0.2 dex. The probability of placing an artificial object at 
a particular location within the image has been weighted by 
the intensity of the smooth component of FUV emission. In 
each case we applied our processing pipeline to the modified 
image. A recovery was considered successful if a region was 
identified within 4 pixels ('set 2') or 2 pixels ('set 1') of the 
original placement with measured flux within 0.2 dex of the 
input 

The 50 per cent completeness thresholds plotted within 
Figure[3]is the flux at which the successful recovery percent- 
age in the Monte Carlo tests is 50 per cent. This value is 
extrapolated between the Ha luminosity bins. We find the 
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Figure 5. Ha luminosity histograms for objects in a given FUV luminosity bin. Histogram colours as in Figure[4l Data histogram i 
only includes Ha selected objects from the clean subsample. As before, data are stacked from the sample of eight galaxies. 



detection threshold within the galaxy to be typically an or- 
der of magnitude higher than for objects placed outside the 
galactic disc. We also find that changing the desired flux 
recovery accuracy to 0.1 dex does not typically increase the 
completeness threshold by more than 0.2 dex. 

An important issue related to crowding is the struc- 
tural complexity of the star formation regions. As a rule, 
the automated object identification algorithms tend to split 
regions with complex morphologies into constituent sub- 
structures. This issue tends to be even more severe in some 
of the oth er investigated object identification software, e.g. 
HIIPHOT (|Thilker. Braun. fc Walterbosll2000l '). While SEx- 
TRACTOR is less prone to splitting the complexes into sub- 
components, the issue is still present in our catalogues, al- 
though to a lesser extent. Bear in mind, this is similar to the 
approach adopted in th e manually compiled catalo gues of 
Hll regions presented bv lHodge fc Kennicutj (| 19831 ), where 
most of the Ha peaks are identified as individual regions. 

Figure [6] shows a field in the south-western part of 
NGC 300. Marked objects fall above the selection line il- 
lustrated in Figure [21 Objects from the Ha selected sample 
are marked in green, the ones from FUV selected sample - in 
magenta. Note that a number of regions have been selected 



independently in both Ha and FUV, with similar spatial 
extent. The majority of objects exceeding the CMDMSM 
constraint are expected to be faint FUV sources. Hence we 
do not expect to be able to identify them reliably using FUV 
selection. This is reflected by the dominance of Ha selected 
objects in our candidate sample. 

The large star forming complex in the southern part 
of the field in Figure [6] is a good illustration of the issues 
with region subdivision. We see a degree of substructure 
present both within Ha and FUV data. The components of 
the complex are treated as individual objects. There are also 
spatial offset between the sources of ionizing radiation and 
the Ha emission peaks within the comples. Ha radiation 
has a clumpy shell structure, surrounding a pair of domi- 
nant FUV peaks. The shell morphology is clearly manifested 
in the Ha/FUV map. The large Ha/FUV values of those 
regions can therefore be explained in those cases by the lack 
of spatial coincidence between Ha and FUV emission. This 
reiterates the importance of using a clean control sample in 
the study. 

An additional qualitative verification with Hubble 
Space Telescope (HST) broadband F606W and F555W im- 
ages of the candidate regions was performed where the data 
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Figure 6. A field wittiin thie SW part of NGC 300. Nortli is up, east is left. Left: Ha image witfi tfie smootfi component subtracted, 
inverted colour scale. Middle: FUV image with the smooth component subtracted, inverted colour scale. Right: Ho/FUV ratio map of 
the two images, obtained by convolving the Ho image to FUV resolution and then dividing it through by the FUV image. Green circles 
are Ha selected regions, magenta circles are FUV selected regions. Note that the Ha/FUV map in the rightmost panel is dominated by 
features originating in areas lacking any physical objects and so is not particularity useful as an identification tool. It is useful, however, 
as means of confirmation of our detections. 



was accesible via Aladin software (jBonnarel et al.ll200Cll V 
This check reveals presence of star clusters at the majority 
of candidate locations. In some cases, three colour compos- 
ites were available via the Hubble Legacy Archive. Example 
morphologies of candidate regions are displayed in Figure 

m 

The object in the upper panel has what could be viewed 
as a model morphology, with a faint shell of Ha surround- 
ing a cluster with a luminous blue star present. The objects 
in the bottom panel have been identified in a much more 
crowded area of the galaxy. The northern region in that 
panel presents a clear cluster morphology. Interestingly, it 
contains bright red as well as blue, which could be evidence 
of multiple, spatially coincident star forming events. The 
southern object also contains luminous blue stars, but there 
is a chance of external ionizing flux contamination associ- 
ated with neighbouring populations. 

To recapitulate, the detections within the outskirts of 
large complexes can be unreliable, especially if wispy Ha 
morphology is seen within the candidate aperture. In those 
cases the ionizing radiation originates outside the aperture 
containing the Hq source. This leads to an artificial boost 
in the observed Ha/FUV ratio. Faint Ha peaks identified 
in relative isolation are considerably more reliable, with co- 
incident star clusters revealed in HST images. 



5.3 Dust attenuation consideration 

In this section we justify the necessity for correcting for 
dust attenuation. We also discuss the influence the choice 
of dust correction recipe has on the results of our study and 
the error balance. 

First, let is consider the effe ct of the chosen correction 
method of ICalzetti et al.l (|2007l ). An alternative approach 
using the relationship between the FUV a nd MIPS 24um lu - 
minosities is also possible, as described bv lHao et al.l l|201lh . 



Finally, it is possible to determine the attenuation values in- 
dependently of each other, which means direction of the at- 
tenuation vector no longer remains fixed by the extinction 
curve. Note that this approach presumes to allow for the 
variation in the relative geometry of dust, stars and ionized 
gas. 

Figure |S] illustrates the results of these three different 
approaches to the derivation of attenuation values. The data 
points are plotted with error bars omitted for clarity. Me- 
dian error bars are indicated in each of the plots. The points 
that are lower limits in Ha/FUV ratios are also upper lim- 
its in I/Fuv • Top left panel shows the uncorrected fluxes 
for ref erence. Top rightpane l illustrates the results of adop- 
tion of Calzotti ct al. (2007) recipe , whil e bottom left panel 
shows the application of Ha o et al.l l|201ll ) relation. The bot- 
tom right panel p resents a comb i ned m ethod where Ana is 
derived from the ICalzetti et al. (12007) pr escription, while 
Apuv is derived from the lHao et al.l I 201ll ) relationship. 



Prior to the application of any corrections, a significant 
number of regions are present at Ha/FUV values exceed- 
ing the theoretical predictions. We note a group of objects 
within the locus of the diagnostic branch of PSS models. Ac- 
counting for dust attenuation effects is necessary in order 
to ensure the physicality of presented populations. 

A striking feature of the plots is the convergence of data 
points in the combined method (bottom right panel) to a 
constant flux ratio value of log(LHa/iFuv) ~ —2.0. This is 
a direct result of the a ssumptions within th e der ivations of 
the tw o relations. Both lCalzetti et al.l (|2007l ) and lHao et al.l 
(|201l|) have calibrated their empirical relations against in- 
tegrated properties of the galaxies. Therefore, the combined 
method recovers the intrinsic Ha/FUV ratio of the assumed 
underlying stellar population. Note also that the reason why 
this approach converges to a value lower than the median of 
our high-mass models is because we are only investigating 



© 2013 RAS, MNRAS 000, [D-?? 



Ha to FUV ratios in star forming regions of nearby spirals 13 




Figure 7. HST view of two fields in NGC 2403, colours are 
B,V,I+Ha . The green circles are three objects exceeding the 
CMDMSM-predicted Ho/FUV flux ratio values. The image was 
taken from the Hubble Legacy Archive. North is up and east is 
left. The top panel is approximately 21 by 12 arcsec, the bottom 
panel 42 by 22 arcesec. Images were produced using the Aladin 
software. 



populations less than 8 Myr old. Including the older objects 
would reduce this value. 

With the combine d method exclude d, we need to inves- 
tigate the behaviour of iHao et al.l (|201l|) and lCalzetti et al.l 
(120071 ) relations when applied to regions with intrinsic 
Ha/FUV ratios deviating from the galactic average. The 
fundamental issue is the dependence of the IR-derived dust 
attenuation corrections on the intrinsic Ha/FUV ratio of 
the objects. For any region with a stellar population too 
old to produce significant ionizing flux (i.e. older than ap- 
proxim ately 10 Myr in the classical view) the lCalzetti et al.l 
([20071) method will overestimate the effect of dust attenu- 
ation. This is because the age-induced intrinsic lack of Ha 
flux will be wrongly attributed to dust absorption. This ef- 
fect is the origin of the high Apuvipuv and low Ha/FUV 
value regions in the top right panel of Figure [S] 

The question becomes more complex in the case of clus- 
ters with high Ha/FUV flux ratios. Those objects normally 
do not populate the IMF fully so that they are intrinsically 
deflcient in the FUV-generating B-type stars. The FUV ra- 
diation field generated by those stars is normally the pri- 
mary source of heating of the dust within the cluster site. 
Hence, any of the available calibrations becomes problem- 
atic. The changes in the thermal balance of dust due to 
miss ing B stars are expected to cause the ICalzetti et al.l 
(|2007i) calibration to underpredict attenuation values. This 
could lead to false positives within the diagnostic regime. 

An additional layer of complexity is added by the ef- 
fects of geometry. The galaxy- wide calibrations average over 



large numbers of lines of sight within the disc, making sim- 
ple geometry models more applicable. Here we are consid- 
ering individual objects with Ha/FUV balance dominated 
by only a few stars. The relative distributions of the dust 
grains, the ionized gas and the stars themselves become of 
a key importance. Exploring these issues would require a 
complex suite of radiative tra nsfer models (e.g. w ith the 
Cloudy code as described bv lFerland et al.l [19981 ) . which 
remains beyond the scope of this study. 

A variation in geometry corresponds to a change in 
the effective attenuation curve, modyfing the proportion- 
ality between Ana and Apuv- We are fully aware of the 
associated caveats. As demonstrated earlier, however, inde- 
pendent attenuation corrections only recover their implicit 
assumptions. We have no feasible means of infering any- 
thing about the line of sight towards the individual clusters. 
Therefore, we apply a straightforward extinction recipe to 
achieve a best average performance across our samp le. 

T o avoid this issue, we adopt the recipe of Hao et al.l 
(|201ll) . Although the use of an Ha based recipe seems in- 
tuitive given that we are using Ha as reference wavelength 
for region identification, FUV-based calibration has number 
of advantages. First, the question of dust thermal balance 
should be better adressed by the method tied directly to the 
heating radiation field. Second, it addresses the problem of 
missing FUV flux due to the use of Ha based apertures. 
The missing flux will increase the attenuation value. This 
will then remove the region from the diagnostic locus, re- 
ducing the chance of false positives. 

We have also te sted the u se of an alternative FUV 
based calibration of iLiu et al.l (|201ll ) (equatio n 5). This 
recipe follows the same functional form as that of lHao et al.l 
(|201ll) . with the coefficient on the MIPS term changed from 
3.89 to 6.0. It was derived for younger stellar populations. 
We found that the higher attenuation values do not alter 
the location of most data points significantly. The most sig- 
nificant effect is the reduction of the number of regions pop- 
ulating the diagnostic regime. We also observe the highest 
Lfvv data points forming a trend aligning with the direc- 
tion of the attenuation vector, which is an indication of 
excessive attenuations. 



5.4 Discussion of uncertainties 

The error bars presented in Figure [3] indicate the uncer- 
tainties associated with photometric measurements and the 
uncertainty in distance affecting the conversion from fluxes 
to luminosities. Other s ources of uncertainty a re the flux 
calibrations themselves. iKennicutt et al.l (|2008l ) quote the 
median uncertainties in their Ha flux calib rations as 12 per 
cent (^ 0.05 dex). GALEX calibrations (JMorrissev et al.l 
[20071) are expected to be accurate to within 1 per cent 
(~ 0.004 dex), which is negligible when compared to the 
Poisson uncertainties and possible systematic effects. 

There are two potential sources of systematic errors in 
this study - the dust attenuation corrections and the pos- 
sibility of inadequate selection of apertures. In the case of 
dust, we have selected the methodology so as to minimise 
the chance of systematics causing false positives within the 
diagnostic region. We appreciate that the geometry vari- 
ations within the regions can plausibly alter the derived 
attenuation values by more than a magnitude. 
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Figure 8. Ha/FUV flux ratio plots of Ha selected region population in NGC 5236 plotted against FUV luminosity with different 
appro aches to dust attenuation corrections. Top left: no correcti on applied. Top r ight: both Apuv a-^d Ajja derived from lCalzetti et al.l 
1120071 ) relation. Bottom left: both Apuv a-^d Ajja derived from lHao et al.l 1120111) relation. Bottom right: Apuv and Ajja derived inde- 
pendently. Note that independent approach converges to a constant Ha/FUV ratio. Contours are the same as in Figure[3] Black crosses 
represent individual data points with error bars ommitted for clarity. Blue arrows represent limits on the ratios, where applicable. Lower 
limits on a flux ratio correspond to upper limits on Lpuv ■ Black arrows represent a dust attenuation correction vector corresponding 
to Ajj„ = 1 mag. The median errors are indicated in top right corner of each figure. 



The influence of aperture selection is also hard to quan- 
tify. Spatial offsets bet ween FUV and Ha p eaks have been 
observed before (e.g. ICalzetti et al.l |2005|). Crowding in- 
creases the challenge of assigning the correct fluxes to indi- 
vidual regions. The visually selected clean sample has been 
created to address these issues and has shown results con- 
sistent with the results from the full data set. Validation of 
candidate regions against HST photometry has also ensured 
the observed diagnostic subpopulation is not merely an ef- 
fect of a systematic error. Furthermore, any objects where 
only flux limits are available have been excluded from the 
histograms in Figures U and [S] Nevertheless, we acknowledge 
that the systematics originating from the source matching 
and background subtraction can modify the fluxes by fac- 
tors up to 0.5 dex in the most extreme cases. 



5.5 Statistical significance of the results 

Statistical treatment of this type of study is made complex 
by the non-trivial variation of completeness limit with po- 
sition, as well as the intrinsic systematics due to the dust 
attenuation correction and our ignorance of any flux leak- 



age. Given our simplistic galaxy- averaged treatment of the 
completeness effects and the lack of model free parameters, 
we do not require our models to return a nominal high prob- 
ability in a single model statistical test. Instead, we focus 
on the relative likelihoods of the models. 

The most general approach in this case is to perform a 
simple Poisson-type test in the regime of our diagnostic high 
Hq/FUV ratio wing. We take the lowest Lpuv bin of Figure 
|4]with the total population of 46 objects and split it into two 
subgroups according to the dust-corrected Ha luminosity. 
We observe 26 objects with logio{LHa/erg s~^) ^ 36.8 and 
20 objects with logio{LHa/eig s~^) > 36.8. When consid- 
ering the clean region sample, we observe 22 objects within 
the lowest Lpuv bin. those values are 13 and 9 respectively. 
We can see that the proportion of populations within these 
bins remains in excellent agreement. The threshold Lhc 
value has been selected to be close to the maximum value 
allowed for the CMDMSM models. 

We estimate the number counts predicted by our mod- 
els by multiplying the normalized model distributions by 
the total observed number of points. However, we wish to 
stress that the model distributions are very sensitive to the 
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accuracy of the approximate treatment of the completeness 
hmits of our sample. We note that the completeness cor- 
rections have been derived in a galaxy- averaged manner for 
dust-uncorrected region properties. Therefore, we do not ex- 
pect either of the region populations to produce a nominal 
high confidence fit to the observed data sample in terms of 
a strict hypothesis testing approach. Despite quoting the 
probability values, we are most interested in the relative 
probabilities between the models. 

PSS model populations predicts Niow,PSS = 37.8 data 
points within the low Lna bin and Nhigh,pss = 8.1 within 
the high Lhu bin. Assuming a Poisson distribution, this 
corresponds to the probability values of Piow,PSS ~ 0.028 
and Phigh,pss ~ 0.00034 respectively. CMDMSM models 
predict Niow,cmdmsm = 45.0 and Nhigh,CMDMSM = 0.9, 
corresponding to Piow,cmdmsm ~ 0.0015 and negligible 
Phigh,CMDMSM (nominal value of 0(10"^*^)). The com- 
bined formal likelihood ratio favours the PSS models by 
more than ten orders of magnitude. While the exact nu- 
merical values are subject to completeness effects, a strong 
preference for the PSS models is made apparent. 



5.6 Notes on individual galaxies 



5.6.1 NGC 300 



The Ha data presented by iLarsen fc Richtlerl l|l999l ) do 
not cover the whole extent of the optically bright disc 
of this nearby spiral. Our calibration has been per- 
formed by matchi n g flux es of corresponding apertures on 
iLarsen fc Richtleij () 19991 ) and wider field, lower resolu- 
tion i mage presented by iHoopes. Walterbos. fc Greenwaltl 
(119961). The scaling was th en set to match the published flux 
of lKennicutt et al.l ([2008) . We acknowledge that this affects 
the quality of the calibration adversely and may contribute 
to the seen offset from the models. This offset could also be 
ex plained by the lower m etallicity of NGC 300, as reported 
bv lBresolin et al.l (|2009l ). The minimum aperture diameter 
used was 15 pixels. 



5.6.2 NGC 2403 

As explored in Section [J] this paper presents a new 
continuum-subtracted Ha image of this galaxy, following 
an astrometry recalibration. A large subsample of regions 
within the diagnostic regime were verified against HST im- 
ages. This galaxy has larger than average uncertainties on 
individual regions' luminosities, arising primarily due to a 
significant u ncertainty in th e dista nce estimate, following 
the study of iFreedman et al.1 (|200ll ). 



5.6.3 NGC 3031 and NGC 5^57 

The available Ha ob servations for these galaxies taken from 
iHoopes et al.l (|200ll ) are of low resolution. This limits our 
ability to detect the regions that would fall within the di- 
agnostic regime. The high dust content within those two 
targets magnifies the imp ortance of an appropr iate attenua- 
tion correction. Using the lCalzetti et al.l (|2007l ) recipe leads 
to a much larger number of regions being identified within 
the diagnostic regime. 



6 DISCUSSION 

6.1 Comparison to the BPASS models 

We have demostrated earlier that the majority of regions 
within our sample are well matched to the model BPASS 
populations in the high luminosity regime. This is not sur- 
prising, since we expect Ha selection to return mostly young 
star forming regions within the modelled age range of 1-8 
Myr. The regions at these luminosities have higher masses 
(Meci ^ 10^'® Mq) and so we expect them to fully popu- 
late the IMF. Stochastic effects become unimportan t in this 
mass regime, as demonstrated previously by Ko da et al.l 
(l2012l) . 

We also observe a small number of objects with 
Ha/FUV ratios lower than our predictions. These are 
seen clearly e.g. in figure [S] The objects in question 
can be explained by two effects. First, we are not in- 
corporating the effects of flux leakage into our mod- 
els. This is especially significant for the Ha luminosi- 
ties of our regions - in practice, the fraction of ioniz- 
ing radiation escaping from the Stromgren sphere volume 
into the surrounding space may exceed 50 per cen t (e.g . 
Oev fc KennicuttI [l997l: IZurita, Rozas, fc BeckmanI l200ll: 



Relaiio. Peimbert. fc BeckmanI I2OO2I : lEldridge fc Reland 
201H ). The leakage of ionizing radiation could explain 



the lowered Ha/FUV ratios of our objects (see also 
iRelaiio et aL|[2012 ). 

The second explanation is age of the regions. Given 
the 1-8 Myrs age range of our models, we expect clusters 
containing stars older than this upper age limit to have 
lower Ha/FUV ratio values. Star clusters with ages below 
100 Myr are expected to retain significant Lpuv output 
after the massive O stars producing the majority of ionizing 
flux have completed their lifecycles. 

The detection of this subpopulation is to be expected 
within our setup. Selection based on Ha luminosity limits 
the number of regions detected at low Ha/FUV ratio values. 
Adopting the FUV selection reveals a much larger number 
of objects within this locus. The full description of the FUV 
selected sample lies beyond the scope of this paper, since 
it offers no insight into the problem of the mmax — Med 
relation. This topic is going to be investigated in a future 
publication (Hermanowicz fc Kennicutt, in prep.). 



6.1.1 Nature of high Ha/FUV ratio objects 

In this study we detect a subsample of star forming regions 
within the nearby galaxies that appear to be inconsistent 
with the CMDMSM models drawn using a truncated stellar 
IMF. Our models suggest those objects are young stellar 
clusters of low masses (Med < 10^'^Mq), which nevertheless 
contain a massive WR or O type star with initial mass in 
excess of 2OM0 . 

Here we explore the alternative explanations for the 
presence of the high Ha/FUV ratio population. First, these 
objects may all be claimed to be results of different sys- 
tematics. The misalignment of Ha and FUV peaks may 
indeed cause elevated Ha/FUV flux ratios in some areas, 
as we have discussed before. Nevertheless, the clean region 
sample, compiled precisely in order to exclude systemat- 
ics of this nature, displays the same departures from the 
CMDMSM predictions. 
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Alternatively, Hq brightness peaks within the diffuse 
ISM could present as objects with high Ha/FUV values. 
However, HST photometry reveals the isolated Ha peaks 
to be typically associated with faint stellar clusters. This 
proves that the 5cr selection cutoff is normally sufficient to 
eliminate Ha brigthness ffuctuation artefacts from our sam- 
ple. 

A number of published studies report existence of mas- 
sive stars within very low mass clusters. Lamb et al. (2010) 
identified a number of candidate associations within the 
Magellanic Clouds dominated by a single O-type stars. 
Their estimates for the masses of these clusters and their 
dominant stars seem to be incompatible with the presence 
of a stron g cluster mass depen dent IMF truncation. Fur- 
thermore, iBressert et al.l (J2012f ) recently presented a VLT- 
FLAMES Tarantula Survey results indicative of O-type 
stars forming in isolation or within very low mass associ- 
ations. 



6.1.2 Potential for runaway star contamination 

Throughout the study we implicitly assume that any stars 
within the apertures have formed m situ. However, our sam- 
ple could be contaminated by high velocity runaway stars 
ejected from massive clusters. Hll regions associated with 
such runaways would have very little corresponding FUV 
emission and hence would naturally display high Ha/FUV 
flux ratios. We have no information on the kinematics of 
these objects, so we cannot exclude a runaway scenario di- 
rectly. 

Similar concerns ha ve been raised rega rding the stellar 
associations identified bv lLamb et al.l (|20ld ). The stellar ve- 
locity structure of those objects is unknown, it is therefore 
possible that their most massive stars are in fact runaways 
from larger clusters. The neighbouring low-mass stars may 
be merely at a coincident location and no t a part of the same 
star for mation event. A recent study bv iGvaramadze et al.l 
l|2012bl ) presented a bow shock detection around one of the 
stars in that sample as evidence of its runaway nature. The 
same study presents evidence s upporting a runaway nature 
of the remaining stars i n the Lamb et al. (|2010i l sample, 
based on spectr oscopv of [ Massey et al.l (|2009l ) . In the case 
of the sample of IBressert et al.l (2013), the radial velocities 
of the target stars have been measured and the probabil- 
ity of the observed objects being line-of-sight runaways has 
been greatly reduced. However, they could still be runaways 
within the plane of the sky or bound stars that are close to 
the apocentre of a highly eliptical orbit. 

The velocity distribution of runaway stars is 
somewhat uncertain. Runaway status is often diffi- 
cult to verify observationally (see review by IStond 

liggil). espec ially for lower ejection velocities . Works 

of iLeonardI (|l99ll ) and ISilva fc Napiwotzkil (|201ll ) 
find OB runaway s tar velocities as large as 700 km 



Howe ver. iHoogerwerf. de Bruiine. fc de Zeeuwl 



liool) and iTetzlaff. Neuhauser. fc Hohld (|201ll ) find 
that the majority of stars have runaway velocities 
slower than 100 k m s~^. Si mulation s of t he run away 



lower velocities. Therefore, the occurance of massive high 
velocity runaways remains possible, but is nevertheless 
rare; 

[Anderson et al.l (|2012l ) in a study of similar resolution 
to ours, found that 80 to 90 per cent of type Ib/c SNe occur 
in regions associated with Ha emission in their host galax- 
ies. The progenitors of those SNe are belie ved to be WR 
stars(e.g. IPodsiadlowski. Joss, fc Hsulll992r ). Furthermore, 
the distribution of the number of SNe closely follows the 
distribution of star-formation. The remaining 10 to 20 per 
cent occur in areas with no significant Ha emission. The 
natural reduction in Ha luminosity and cluster mass dis- 
tribution implies that some SNe will occur in regions with 
Ha emission below the detection limit. Therefore this study 
suggests that very few of the WR stars are runaways that 
travel a large distance from their birthplace. 

Finally, it is uncertain whether any potential WR or O 
type runaways will produce observable Ha emission. While 
every O or WR runaway produces ionizing photons, the 
density of the surrounding gas determines whether an H ll 
region is observed. Unless the local density is high, the H ll 
region might be extremely large and therefore of too low a 
surface brightness to be detected in our sample. If sources 
we have detected are associated with runaway stars, they 
must but be in a region of high ISM density so that the Ha 
emission is constrained within a compact region less than 
^ 70pc in diameter. This corresponds to the ISM number 



density of the order of ~ 1 cm for 09 stars and 



10 



star population by Fuiii fc Po rtegics ZwartI (|201l|) and 
lEldridge. Langer. fc Touti l|201ll ) also show that massive 
runaways are typically ejected with velocities significantly 
lower than 100 km s~^ and that more massive stars attain 



cm~ for 03 stars in the Str omgren sphere approximation 
IjOsterbrock fc Ferlandll2006l ). Th ose values are with in the 
regime of giant molecular clouds (jHeiner et al.ll2008l ) even 
prior to accounting for geometry effects and the ambient ra- 
dia tion field. Objects like this are known (e .g. C, Ophiuchii, 
see IGvaramadze. Langer. fc MackevI l2012al . and the refer- 
ences therein), but expect ed to be rare. Fur t hermor e, as- 
suming the proposi t ion of iGvaramadze et al.l (|2012b) that 
all the JLamb et al.l l|2010r) O type stars are runaways, we 
note that only two of the stars have surrounding Ha emis- 
sion that could be discernible at the distances of the galax- 
ies analysed in the paper (b ased on the Magel lanic Cloud 
Emission Line Survey images ISmith et al.ll2000r ). 

We also note that under the scale- free star formation 
paradigm proposed bv iKruiissenI (|2012r ). most star forma- 
tion in spiral galaxies is expected to be in unbound asso- 
ciations. Those slowly dissolving structures do not undergo 
dynamical relaxation and are not expected to produce high 
velocity runaways via dynamical interactions or violent gas 
expulsion. 

We expect objects populating our high Ha/FUV ra- 
tio branch to be of similar nat ure to the stel l ar ass ocia- 
tions identified in the studies of IBressert et al.l (|2012l ) and 
iLamb et al.l (|2010f ). Massive runaway stars are a poten- 
tial contaminant, but considering all the factors discussed 
above, we do not believe that they account for all the re- 
gions in the diagnostic wing. 



6.2 Stochasticity and the truncation of the IMF 

This work is complementary to a immber of previous studies 
that attempted to verify the presence of a cluste r mass de- 
pendent truncation in the stellar IMF. iLee et al.l l|2009r ) at- 
tempted to observe the Ha/FUV ratio downturn proposed 
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bv lPflamni-Altenburg et al.l l|2009f) in the integrated fluxes 
of galaxies in the local Universe. They found a small de- 
chne in the Ha/FUV ratio at low Lfuv values, in agree- 
ment with the semi-analytic IGIMF predictions. However, 
they found the s ame d ecline to be consistent with models of 
Tremonti et al.l (^2007^ ■ which assumed stochastic sampling. 



Fumagalli et al.l (|201l|) extended this analysis with the use { 



of their SLUG stochastic population synthesis code. They 
have demonstrated that the truncated IMF models signifi- 
cantly underpredict the integrated Ha luminosi ties of galax- 
i es in the lowest FUV luminosi ty bins. Finally, IWeisz et al.l 
(J2012l ) extended the analysis of lLee et al.l (J2009r ) by includ- 
ing the estimated masses of the individual galaxies into their 
analysis, finding that implementation of truncated IMF sig- 
nificantly underpredicts the Ha/FUV ratios of the low-mass 
end of their sample. They offer time variations in the galac- 
tic star formation rate as an alternative explanation of the 
observed sample behaviour. We note their study does not 
account for the effects of stochasticity, though. 

A few studies similar to ou r proj ect have been per- 
formed previously. ICorbelli et al.l (|2009l ) presented a study 
of star clusters in M33, focusing on the relation between 
the Ha and bolometric luminosities. Their study revealed 
existence of clusters with Ha fluxes in excess of what could 
be explained by models with a cluster mass dependent IMF 
truncation imposed. Our study presents evidence of pres- 
ence of analogous objects in a number of more distant spiral 
galaxies. We also base our conclusions on an updated set of 
models, which include the effects of stellar evolution. 

ICalzetti etHI (|2010l ) used the HST to search for such 
objects within NGC 5194. Their investigation revealed a 
population of clusters (with 3a detections in Ha) consis- 
tent with fully stochastically sampled IMF, with no evi- 
dence supporting a presence of a cluster mass dependent 
IMF truncation. The results become less conclusive when 
the full object sample is taken into account. The study is 
based on a HST BVI-selected cluster sample des cribed in 
detafl bv lChandar et all l|201ll) . lKoda et all (|2012l ) have re- 
cently presented a study of star formation in the XUV disc 
of M83 (NGC 5236), demonstrating their results to be con- 
sistent with stochastically sampled Salpeter IMF. We find 
our conclusion s agree with the results of both those studies. 
ICalzetti et al.l (|2010r ) focused on comparing the Ha lumi- 
nosity of the clusters to masses extrapolated from U,B,V,I 
photometry. Our approach incorporates the directly ob- 
servable FUV luminosities, he nce directly addressing the 
IPfiamm-Altenburg et al. (|2009|) predi ctions in the regime of 
individual clusters. iKoda et al.l (|2012l ') focused on the num- 
ber counts of Ha - and FUV-bright clusters with a sim- 
ple ageing relation explaining their relative number counts. 
Their investigation of Ha/FUV flux ratios was limited to 
the integrat ed properties of t he XUV disc, being parallel to 
the study of iLee et all (|2009l ). 

In all of those cases, our study complements the results 
by expanding the investigated sample of regions and allow- 
ing for attenuation variation between the individual objects. 
We also base our conclusions on an updated set of models, 
which include the effects of binary interactions on stellar 
evolution. Our study flnds apparent disagreement with the 
truncated IMF model predictions within the regime of in- 
dividual cluster s ' prop erties , matching the conclusions of 
iFumagaUi et all (|201lf '). and lEldridgl l|2012[ ). However, the 



observational constraints from the dataset presented in this 
paper are not sufficient to exclude weaker manifestations of 
the IMF truncation. 

Note that our data are not sensitive to differentiate be - 
tween the random sampling preferred by Elmegreenl 1I2OO6II 
and t he sorted sampling preferred by IWeidner fc Kroupal 
20061) . We note those two methods predict very similar 
property ranges for individual clusters. Wc also wish to 
stress that the work of PA07 and Pflamm-Altcnbur g et al.l 
l|2009r ) are based on semi-analytic approximations and 
therefore unsuitable for stellar population synthesis using 
discrete stellar models. The PA07 relationship was derived 
to represent the typical rrimax at a given M^d and stochas- 
tic departures from this relationship are expected, both to- 
wards higher and lower stell ar masses. 

In line with the work of ICervifio et al.l (|2003l ). we con- 
clude that the use of a fully stochastic stellar population 
synthesis code is critical in any study of small mass clus- 
ters and low SFR galaxies o r environments. Codes like the 
widely used Starburst99 (JLeitherer et al.l 1 19991 ') are in- 
herently unsuitable for the investigation of the effects of 
stochasticity. Fixed implementation of the IMF means that 
these codes are unable to represent the diversity of prop- 
erties displayed by low mass clusters and young galaxies 
with low SFRs. We suggest that stochastic codes are used 
in modelling of stellar populations in such cases. 



6.3 Caveats 

There is a number of effects that are not accounted for in 
this study. This section is a short overview of their relative 
importance. 

Our analysis makes no provisions for the dynamics of 
the systems in questions. One such a scenario could be a 
formation of a close binary in the centre of a low mass 
cluster consistent with the IMF truncation, ejection of a 
signiflcant portion of the cluster mass (probably mainly via 
ejection of the other massive stars) and then formation of 
a massive star with m > mmax via a coalescent binary sce- 
nario. Whilst we allow for effects such as common envelope 
evolution and mergers within the scope of our binary mod- 
els, multi-star dynamical interactions are not included in 
the models. In practice, the stellar merger rate could be 
signiflcantly boosted due to dynamical interactions within 
the cluster. Hence, we would see more CMDMSM-consistent 
clusters populating the diagnostic high Ha/FUV ratio area 
in Figure (21 This would require an increase of the merger 
rate by two order of magnitude. 

However, the recen t numerical simulations presented 
by lOh fc Kroupal l|2012t ) suggest that the effects of dynam- 
ical evolution on the rrimax vs. M^d relation is negligible. 
They flnd that only a small fraction of low mass embed- 
ded clusters will have the most massive star's mass signif- 
icantly altered by stellar collisions. The quoted study only 
considered clusters formed under IGIMF paradigm. In view 
of these results, we do not expect the dynamical evolution 
to be responsible for the presence of high Ha/FUV ratio 
population in our data. 

A related concern is the potential effects of violent gas 
expulsion on the evolution of the observed clusters. Numer- 
ical s tudies (e.g. iBaumgardt fc Kroupal l2007l : ISmith et al.l 
I2OI3I , and the references therein) suggest that it is not un- 



© 2013 RAS, MNRAS 000, [T]-?? 



18 M.T. Hermanowicz, R.C. Kennicutt and J.J. Eldridge 



common for the clusters to lose even 90% of their associ- 
ated gas mass and as a result a significant population of 
stars formed in the original star formation event can be- 
come unbound. Smit h et al.l (|2013l ) point out that the asso- 
ciated dynamical evolution of young clusters is a stochastic 
process highly sensitive to the initial conditions. However, 
under such a scenario we would expect the majority of the 
unbound stars not to have travelled far enough from our 
clusters so as to fall outside the apertures used in this study. 

The BPASS models do not include the effects of as 
rotation in stars, X-ray binary systems, or wind-wind col- 
lisions within the binary systems. The inclusion of any of 
these would lead to an increase in the ionizing radiatio n rel- 
ative t o FUV. In fact, recent models bv iEldridge fc Stanwavl 
()2012l ) indicate that incorporation of rotation into massive 
binary star models lead to a 30 per cent increase in the 
ionizing flux. 

We are also assuming that the whole extent of the FUV 
radiation is originating solely from the underlying stellar 
population, hence ignoring any potential contributions from 
the ionized gas itself. Even though ignoring this contribu- 
tion would arti ficially boost the Hq/ FUV values of our ob- 
jects, work bv lDanforth et al.l l|2003l ) indicates that even in 
the presence of shocked gas within a supernova remnant, the 
UV spectral energy distribution (SED) is clearly dominated 
by the stellar continuum. 

In our analysis we do not consider the variations in 
the metallicity. We base our analysis on the models of So- 
lar metallicity, where i n reality the metalliciti es vary from 
one galaxy to another (JMoustakas et al.ll201(J l. Additional 
radial trends have also been observed within the discs of 
spiral galaxies (e.g. iBresolin et al.ll2009l . and the references 
therein). Due to the issues with crowding, we rarely find any 
high Ha/FUV ratio regions close to the centres of studied 
galaxies. However, where the crowding becomes less of an 
issue, we do not observe a clear preference for identification 
of high Hq/FUV ratio regions at high galactocentric radii. 
The image of NGC 300 used in our study only covers the 
inner part of the galactic disk, and yet objects inconsistent 
with CMDMSM are identified there. 



7 SUMMARY AND CONCLUSIONS 

We detect a small subpopulation of regions with high 
Hq/FUV flux ratios. Comparison with model populations 
suggests those regions to be consistent with the results of 
stochastic s ampli ng of an universal IMF, as suggested by 
lElmegreenI (l200a ). We note that sorted sampling and PSS 
model populations are not easily distinguishable in the type 
of study presented here. 

We find that any study of those objects within the 
galactic discs is sensitiv e to the ad o pted dust attenuation. 
An FUV-based recipe of lHao et al.l (|201ll ) was found to be 
most unbiased when combined with the effects of stochas- 
tically populated mass functions within the clusters. The 
high Ha/FUV ratio population persists after an applica- 
tion of dust attenuation corrections. 

The candidate region population is tested directly 
against the predictions of truncated IMF cluster models. We 
identify a subsample of objects whose properties seem to be 
inconsistent with a cluster mass dependent IMF truncation 
in the form described in Section 14.11 This view is consis- 



tent with a number of previous studies (ICorbelli et al.ll2009l : 

;|k " 



ICalzetti et alll2010l : IWeisz et al.ll2012l : iKoda et al.ll2012l ). 

We emphasise the significance of stochastic effects in 
the studies of low SFR environments and caution against 
the application of semi-analytic recipes and classical stellar 
population synthesis code in this regime. 

One of the factors limiting the scope of this work 
is the low resolution of GALEX data available. Analo- 
gous studies will become more definitive with the future 
launch of World Sp ace Observatory-UV planned for 2016 
l|Shustovet al.ll201ll ). 
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APPENDIX A: RESULTS FOR CALZETTI 
ATTENUATION CORRECTION 

In this appendix we present alternative versions of the Fig- 
ures [3], [4] andO which have been created with the applica- 
tion of the Calzetti ct al. ( 200' g) dus t attenuation relation, 
as opposed to the iHao et al.l (201l|) relation used in the 



main analysis. Note the higher scatter of objects around the 
model predictions, especially in the high luminosity regime. 
A larger number of objects is also observed in the diagnos- 
tic wing, making the constraint a gainst the t runca ted IMF 
stronger. This is because the lCalzetti et al.l (120071) attenu- 
ation relation is not as heavily biased against the regions 
with boosted Hq/FUV ratio. 
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Figure Al. As in Figure [S], but witli ICalzetti et al.1 l|20071 ) attenuation recipe applied in tlie place of iHao et alJ l|201U ). Note the 
presence of larger scatter at high luminosities. 
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Figure A2. As in Figure |4], but witli lCalzetti et alj 1120071) attenuation recipe applied in the place of lHao et alj 1I2OIII ). 
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Figure A3. As in Figure[5], but witti lCalzetti et alj ll2007h attenuation recipe applied in tlie place of lHao et alj 11201 Ih . 
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